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Critical opalescence is a characteristic experimental signature of a second order phase transition 
O ' in solid state physics. A new, experimentally accessible measure of opacity and of attenuation 

' length in heavy ion reactions is suggested, as a combination of HBT radii and nuclear modifica- 

tion factors. This opacity is maximal when ^/snn, the system size and centrality correspond to 
■"sj" ■ the critical point of QCD. To characterize the phase transition at this critical point, the critical ex- 

■ ponent of the correlation function can be determined by measuring the Levy index of stability of 

\^ , the Bose-Einstein or HBT correlations. The exponent of the correlation length can be determined 

from fits to the multiplicity distribution in various pseudorapidity intervals, also as a function of 
colliding energy, system size, centrality and (chemical) freeze-out temperature. These two critical 
' exponents determine the remaining four critical exponents and the universality class of this sec 

0^ ■ 



ond order phase transition. As a control experiment, the determination of the critical exponent of 
the specific heat capacity is proposed, from event-by-event fluctuation measurements. To measure 
opacity precisely, well calibrated high transverse momentum probes are needed, such as given by 
the excitation function of y-n jet correlation functions. 
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1. Introduction: five important milestones in Au+Au collisions at RHIC 

The RHIC heavy ion program created a paradigm shift in the exploration of the phase diagram 
of strongly interacting, hot and dense hadronic matter: instead of discussing a gas-like, weakly cou- 
pled plasma of quarks and gluons, we now explore the properties of a perfect fluid, where valence 
quarks appear as active degrees of freedom, and we discuss superfluidity at extremely high temper- 
atures on the scale of T ~ 2 x lO'^ K. An important first milestone in this program was a PHENIX 
discovery of suppression of particle production with large transverse momentum in central Au-i-Au 
collisions at ^Jsj^ =130 GeV, which apparently was the first glimpse on a new phenomenon 
never seen before in high energy nuclear collisions. Soon the effect was found at higher energies 
too [0], and in a control d-i-Au experiment all the four RHIC collaborations reported the absence of 
suppression of particle production at high transverse momentum [^, 0, ^ The angular correla- 
tion plot of STAR indicated particularly clearly, that the most energetic, near-side jets have similar 
angular distributions in p-i-p, d-i-Au and Au-i-Au collisions, while the momentum-compensating 
back-to-back jets are seen only in p-i-p and d-i-Au collisions, but are absent in central Au-i-Au colli- 
sions. This absence of the back-to-back or away-side peak but the presence of the near-side peak 
indicates that jet suppression in Au-i-Au collisions is not an initial state phenomenon, but a final 
state effect. These results convincingly demonstrated that a new form of matter is seen in central 
Au-i-Au collisions, that is so opaque that it attenuates even jets of particles with high transverse 
momentum, and can be considered the second milestone of the RHIC program. The third mile- 
stone has apparently been the publication of "White Papers" by all the four RHIC collaborations, 
that summarized the conclusions from the first three years of data taking at RHIC, and pointed out 
that the matter created in central Au-i-Au collisions does not behave like a gas of weakly interacting 
quarks and gluons, instead it flows like a perfect fluid of strongly interacting quarks and gluons 
[0, [s], IC]. This consensus result later on became known as the discovery of a perfect liquid at 



RHIC, and it has been selected as the Top Physics Story of the American Institute of Physics in 
2005, the World Year of Physics |[Tl|]. Among the many additional milestones to come, let me 
highlight two, particularly important discoveries, and consider the discovery of the quark number 
scaling of the elliptic flow as Milestone 4. This discovery can be interpreted in such a way, that 
the flowing matter has quark degrees of freedom. STAR has pointed out an approximate quark 



number scaling of the elliptic flow in terms of the transverse momentum pt [12], while PHENIX 
strengthened this argument significantly by pointing out that the quark number scaling becomes 
perfect at low transverse momentum if instead of pt, the transverse kinetic energy KEj = mj —m 
is used [|l^] (where the transverse mass is introduced as m, = ^Jm?- + pj). At mid-rapidity, this 
transverse kinetic energy is a reasonably good approximation of the (dimensionless and in general, 
rapidity dependent) scaling variable w, found in first in exact solutions of perfect fluid fireball hy- 



drodynamics 1 14]. When plotted in terms of the universal scaling variable w, all presently known 



elliptic flow data at RHIC collapsed to the theoretically predicted, universal scaling curve []15|], con- 
firming the perfect fluid picture. As the relationship between the scaling variable w and the elliptic 
flow V2 is approximately linear for V2 <C 1, this universal scaling law is approximately invariant for 
a rescaling by the number of quarks. By now, this transverse kinetic energy scaling of with the 
number of valence quarks has been extended to even charm and strange mesons and baryons. The 
scaling of the elliptic flow with the number of valence quarks indicates experimentally, that in this 
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perfect fluid valence quarks are relevant, active degrees of freedom. Such an active role of quarks 
has been expected, but it has not been demonstrated experimentally before [12, 13]. This is why I 
consider this observation as the fourth important milestone at RHIC. 

As a next step, it was natural to ask [16], how perfect is this perfect fluid? Can we measure 
its kinematic shear or bulk viscosity? What about its heat conductivity? The first answer to this 
question was surprising to many: The measured kinematic shear viscosity [ jl7| ] in Au+Au colli- 
sions at RHIC proved to be significantly below that of the same value for superfluid ^He at the 
A -point [18, [T^. At the same time, the temperatures measured in Au-hAu collision at RHIC were 
the hottest temperatures ever measured in experiments, of the order of 2 x 10^^ K, in sharp contrast 
to the ultra-cold temperatures of superfluid ^He, of the order of 4 K. In this sense, we may discuss 
high temperature superfluidity as the fifth important milestone discovered in Au-i-Au collisions at 
RHIC [|^]. 

These discoveries and important milestones are far from being a complete and exhaustive 
summary of discoveries made in the RHIC proton-proton and heavy ion program. They were 
highlighted here in order to demonstrate, that a new form of matter is created in heavy ion collisions 
at ^Js}^ = 200 GeV energies. This matter is so dense and so opaque, that it absorbs about 80 
% of the jets with high energy; it has temperatures that are higher, than any other temperature 
measured before; it exhibits high temperature superfluidity i.e. flows within errors perfectly; even 
the relatively rare and heavy strange and charm quarks flow with this matter. These properties are 
apparently absent at lower energies, so it is natural to ask, at what bombarding energies, what type 
of collisions and at what kind of centralities do these phenomena show up first, and what are the 
characteristic properties of these transitions? Can we map out the phase diagram of QCD? 

Various theoretical considerations suggested that at high baryon densities there exists a line of 
first order phase transitions on the QCD phase diagram, and this line ends in a critical end point 
(CEP), where the phase transition changes from a first to a second order transition. At even smaller 
baryochemical potentials (or net baryon densities) the transition is expected to become a cross- 
over, i.e. a rapid increase in the number of degrees of freedom but without a clear thermodynamical 
separation of the confined and the deconfined phases. These ideas are illustrated in Fig. |l|. 

Hence it is of fundamental interest to locate the critical end point of the line of first order 
QCD phase transitions: this point acts as a landmark or a sign-post on the QCD phase diagram, 
and separates the line of first order phase transitions from a cross-over region [0]. The next 
section discusses how to locate, identify and characterize such a critical point in high energy particle 
and nuclear physics, based on analogies with similar experiments in physical chemistry and on 
applications of the theory second of order phase transitions. 



2. Critical Opalescence: Signature of a 2nd Order Phase Transition 

Before entering a discussion on how to characterize a critical point, let us first discuss a strategy 
how a critical point is located in statistical physics and in experiments on second order phase 
transitions involving ordinary gases and fluids. This subject will lead us to a new, experimental 
definition of opacity in high energy heavy ion physics. In the next section we discuss how to 
characterize such a critical point by critical exponents, based on the theory of second order phase 
transitions in statistical physics and the adaptation of these results to heavy ion physics. 
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Schematic Phase Diagram of Hot and Dense Strongly Interacting IVIatter 



> 

0) 



Quarl< Gluon Plasma 
(color deconfined) 



Critical Point 

Line o*-^ 

hadron gas 
(color confined) 





nuclear superfluid 



color-flavor 
locked 



-310MeV 



Figure 1: Schematic, theoretically expected phase diagram of QCD. The critical point is shown as the end 
of the line of first order deconfinement phase transitions, based on ref . [ 21 1 . Ref . [|2|] predicted Ie = 1 64 ± 2 
MeV and /ig = 360 ± 40 MeV for the location of the end point of the line of first order QCD deconfinement 
phase transitions. An earlier prediction of ref. [ |23| ] suggested that the critical point is approximately at twice 
as large baryochemical potentials, suggesting Te = 160 ± 3.5 MeV and /X£ = 725 ± 35 MeV. 



Critical opalescence is a characteristic phenomenon that arises in second order phase transi- 
tions, as illustrated in Fig. ^ It is commonly demonstrated in binary fluid mixtures, or at the end 
point of the line of first order liquid-gas phase transitions. As the critical point is approached, the 
sizes of the gas and the Uquid region begin to fluctuate on increasingly large length scales, eventu- 
ally including the size comparable to the wavelength of light. These fluctuations then scatter light 
and cause the normally transparent fluid to appear cloudy or opaque. At the critical point, density 
fluctuations appear on every length scale, in ordinary fluids these scales include the wavelength of 
light (from 700 nm to 400 nm) up to as much as a few mm. The increased opaqueness of the fluid 
is a well visible signal of the vicinity of a critical point. 

In optics, the opacity K is defined as the rate of absorption or extinction, at which the intensity 
/ of a beam of radiation is absorbed or scattered per unit distance along a ray of propagation: 

dl 

^ = -Kl (2.1) 

ax 

The opacity K may depend on the wavelength (or momentum) of the radiation, as well as on the 
density, temperature and composition of the medium. For example, air has nearly zero opacity for 
visible light, and for radio waves, but it is almost completely opaque for gamma and X rays, and 
in most of the infrared part of the spectrum. The solution of the above differential equation for 
constant values of opacity is 

/ = /oexp(— K".)c) = /oexp(— x/A) (2.2) 
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Figure 2: Locating the critical point in physical chemistry by the phenomenon called critical opalescence. 
In a second order phase transition, scattering centers (density fluctuations) appear on all scales, hence the 
matter becomes opaque, as indicated by increased scattering of a penetrating probe, in this case a laser beam. 
The maximum of the opacity corresponds to the critical point. Pictures from ref. [^sj]. 

where A = 1 / k: is the so called attenuation length or penetration depth, defined as the distance 
where a radiation inside a material falls to\/e 36.7% of its incoming intensity. So a characteris- 
tic signature of critical opalescence is a minimum of the attenuation length or a maximum value of 
opacity on many wavelengths where the matter was previously transparent. Measuring opacity as 
a function of pressure and temperature, the critical point can be located. The location of the critical 
point {TcPc) on the {T,p) diagram is an important characteristics of the given material. 

2.1 Opacity in heavy ion collisions 

It is interesting to rewrite the defining relation of optical opacity in a manner that has a straight- 
forward generalization to heavy ion physics: 

/(generated) — /(transmitted) 

K* — — , 

/(generated)Ajc 

where Ax stands for the distance traveled in the medium. 

In heavy ion physics, the nuclear modification factor Raa can be defined as 

^ /(transmitted) /(measured) 
/(generated) /(expected) ' 

where the measured yield is defined as 

while for point-like processes like production of high transverse momentum jets the generated rate 
is the same yield measured in proton-proton reactions scaled up by {Ncou), the number of binary 
initial nucleon-nucleon collisions, given by a Glauber model calculation: 

/(expected) = -^-^. (2.6) 
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This nuclear modification factor Raa includes a mixture of hot and cold nuclear matter effects [p^], 
its most apparent shortcoming is that it measures the change in intensity but does not specify 
the distance over which the change of the yield happened. Similarly, the nuclear modification 
factor can be smaller for a matter which has larger opacity, if this matter is present only in a 
short distance along the line of propagation of the probe, as compared to a less opaque matter 
that has been created in heavy ion colUsion in a larger volume. Opacity as defined in optics, on 
the other hand, compares the rate of change of intensities per unit distance. Thus for a proper 
definition, both the nuclear modification factor and the characteristic length scale of the attenuation 
has to be included. Let us consider, for simplicity, a one dimensional symmetric distribution with a 
characteristic radius of R. Suppose that a pair of back-to-back jets is created inside this distribution 
and the near-side jet can escape after covering a distance r. Obviously, the far-side jet has to 
propagate over a distance of 2R — r. Averaging this over the source distribution with the distribution 
of the production points r the near side jet travels (r) while the far side jet travels 2R — (r) distance 
in the medium on the average. Averaging over both jets the mean distance of propagation becomes 
just R, the characteristic radius of the source. 

Such a characteristic source size is measured from two-particle correlations using various tech- 
niques, named as femtoscopy in general. More specifically, we measure the Bose-Einstein corre- 
lations of like-sign pions and from this measurement in a Gaussian approximation to the source 
density distribution effective source sizes are extracted. These radii are frequently referred to as 
Hanbury Brown - Twiss or HBT radii to honor the two radio astronomers, who invented a similar 
technique for photons to measure the angular diameter of main sequence stars. These HBT radii 



can be (and apparently are [27] ) measured as a function of the angle with respect to the reaction 
plane. For a small attenuation or 1 — Raa ~ £, we may define the measurable opacity as Kg = 
where we divide the relative attenuation of intensity by Rhbt, the characteristic length scale of the 
fireball measured (in principle) along the line of propagation. (The application of Bose-Einstein or 
Hanbury Brown - Twiss correlation techniques are detailed in subsequent sections.) 

It is important to note that in high energy heavy ion reactions in the pt > 4 GeV kinematic 
range we are measuring not tiny absorbtions, but quite significant values, so the thickness of the 
absorber or the fraction of the absorbed intensity cannot be considered small, 1 — Raa ~ 0.8. Hence 



we have to solve eq. ( |2.2[ ) to get a more precise definition, keeping in mind the interpretation of 
the nuclear absorbtion factor as Raa =l/h- Hence for arbitrary thickness and nuclear modification 
factor, the optical opacity should properly be defined as 

,= J^ (2.7) 
Rhbt 

The opacity k defined above may depend in principle on system size, colliding energy, central- 
ity, and angular direction with respect to the reaction plane, through Raa as well as Rhbt- As the 
average HBT radius as well as Raa are both positive, the above quantity is well defined. Using the 
above definition, we can compare the opacity of different media that might be available in differing 
amounts. This definition has several other advantages: for < Raa < 1 the opacity is positive as the 
HBT radius is always positive. In case of no nuclear modification, Raa = 1 and K" = 0, the opacity 
vanishes, indicating that the medium is fully transparent. On the other hand, if Raa > 1> the opacity 
becomes negative, indicating that the matter is not absorbing but in fact emitting extra radiation. 
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However, the main advantage of the above definition is that it defines opacity as a modification 
of spectra per unit length along the direction of the propagation of a penetrating probe. Thus it 
characterizes a local property of the medium, and not an integral property. Local and not integral 
properties are relevant in a search for a phase transition: it may happen that the medium created in 
central Cu + Cu collisions has less absorption or has a larger Raa than that of a mid-central Au +Au 
collision simply because the bulk amount of the medium is larger in mid-central Au +Au collisions. 
Comparing to the HBT radius is necessary to correct opacity for the bulk effect. 

It is even more intuitive to introduce the inverse of opacity, as the attenuation length: 

^ = \ = 'wty ^'-'^ 

This measures the length over which the intensity of radiation is attenuated by a factor of \/e. In 
case of some nuclear absorption, < Raa < 1> and < X < °°. For no nuclear modification and 
Raa = 1> the attenuation length diverges and it becomes negative for Raa > 1, as the medium en- 
hances the radiation. Thus a negative A means that the intensity decreases in the direction opposite 
to that of the propagation. 

So in principle, one has to infer the actual distance that the probe had to cross before departing 
towards the detectors. It is well known that the size of nuclei change on the average as A'/^, which 
in a collision would then be modified to an average Np^art dependence. The freeze-out volume (more 
precisely, the size of the region of homogeneity) of charged pion emission was measured in 200 
GeV Au-i-Au collisions using Bose-Einstein or HBT techniques as a function of centrality [p7|]. 

1/3 

The Gaussian HBT radii were found to scale as /?,■ = po + Pi^parf where po « 0.83 ±0.25 fm and 
pi 0.54 lb 0.05 fm and / = (long, out, side) stands for the three principal directions with respect 
to the beam and the mean transverse momentum of the particle pair. Note that ref. [ p7| ] fitted the 
centrality dependence of each HBT radius components separately, but let me point out that the fit 
parameters in the different directions were found, within one standard deviation, to have the same 
values. Thus their average value was quoted above. Experimentally, 7?o„f « ~ Riong ~ Rhbt, 
due to the effective spherical symmetry of the pion emitting source. Note also that penetrating 
probes (like high p, jets, leptons or gammas) can be radiated, in principle, throughout the time 
evolution of the expanding fireball created in high energy heavy ion collisions. So the distance that 
is covered by these probes in the medium is not a static distance, but a length scale of the expanding 
source that is covered by the jets before they punch through. Few exact results are known for the 
time evolution of such mini-bangs or exploding fireballs of heavy ion reactions, but two points 
are quite remarkable: in a non-relativistic approximation it was found that although each principle 
directions of the expanding fireball diverges to infinity, the HBT radii or lengths of homogeneity 



tend to a time and direction independent constant, so the HBT volume becomes spherical [38]; 



second that experimentally such a spherically symmetric HBT volume is seen in pion-correlation 



data ||37|]; third that the HBT volume - just like the initial volume - is found to scale with Np^lf. 



Based on these three observations we may expect that the distance covered by the probe in the 

1 /3 

medium is indeed proportional to Np^^i, hence scaling (1 — Raa) by the HBT radii that are also 

1/3 

proportional to Np^n-f gives the centrality dependence of the opacity K correctly; furthermore the 
offset parameter po > guarantees that even in the most peripheral collisions the initial volume 
that has to be propagated through by the high pt probes is non-vanishing. Such expansion effects 
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may indeed introduce certain absolute calibration error on the opacity K that is difficult to control. 
However, the relative error on the centrality dependence of the opacities K and in particular its 
minimum or maximum structure is not affected by these expansion effects, as both the initial and 
the final volumes scale proportional to Np^an- 

In principle, both the effective source size Rhbt, and the nuclear modification factor Raa 
may change with respect to the reaction plane even in events with the same colliding system, 
bombarding energy and centrality. Then opacity should be defined by the above ratio evaluated as a 
function of angle with respect to the reaction plane. Although both Raa and Rhbt may change with 
respect to the orientation of the reaction plane, the optical opacity is expected to be independent of 
the selected directions in a fireball that consists of one given type of matter. 

As a test of the proposed method, and for simplicity, let me use here an average HBT radius as 
Rhbt = {Rout +Rside +Riong)/^ and also nuclear modification factors of neutral pions at pt = 4.75 
GeV in PHENIX to evaluate the opacity and the attenuation lengths in ^/s^n = 200 GeV Au+Au 
collisions, as a function of centrality. The results are shown in Table 2A . We find that changing 
centrality from 50 - 60% to 0-5 %, the opacity is increased by substantially, by about a factor of 2 
(a 4 a effect), while the nuclear modification factor decreases by more than a factor of 4 from 0.72 
to 0.19. Thus nearly half of the suppression increase originates from the increased amount of the 
hot and dense, strongly interacting matter created in more central collisions, while the attenuation 
length is also decreased by a factor of 2 in these collisions. 



Centrality 


0-5 % 


20-30 % 


30-40 % 


40-50 % 


50-60 % 


Opacity K (fm^^) 


0.35 ± 0.04 


0.27 ± 0.03 


0.26 ± 0.04 


0.12 ±0.02 


0.15 ±0.05 


Attenuation A (fm) 


2.9 ± 0.3 


3.7 ± 0.4 


3.8 ±0.6 


8.1 ± 1.5 


6.5 ± 2.0 



Table 1: Examples of opacities k and attenuation lengths A in ^/snn ~ 200 GeV Au+Au reactions, evaluated 
from nuclear modification factors measured at pi = 4.75 GeV/c in ref. [ |3^ ] and using HBT radii of ref. p7|], 
averaged over both directions and charge combinations, at the same centrality class as Raa- 



For a precise absolute calibration of the opacity k:, additional information is needed both in 
the soft and in the hard transverse momentum sector. To measure the change of intensities at high 
transverse momentum, we have to know precisely what the initial energy of the jet was, before 
it became attenuated. This is possible to learn, if we trigger on(7±^^7±^or^ + g— >^±7 
events with high pt y-s, and measure the energy of the high transverse momentum photon precisely, 
because y's are interacting only electromagnetically so they are not attenuated in nuclear matter. 
With the help of such processes, the energy and momentum of the momentum balancing hadron 
jet, or the initial intensity can be measured precisely. However, about two orders of magnitude 
more events are necessary for such a study, than for a usual Raa measurement, due to the relatively 
small cross-sections and large 7 background from tt" and other hadronic decays. Similarly, the 
time evolution of the size of the expanding fireball could in principle be measured by soft 7 ± 7 
correlations, as these probes penetrate the hadronic fireball too; but up to now no such studies 
were published in the RHIC program, due to difficulties in background rejection and the need for 
extremely large statistic. We leave these topics as future challenges to the RHIC experiments. 
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3. Critical Exponents and the Universality Class of the Critical Point 

Phase transitions that involve latent heat are called first order phase transitions. During such 
phase transitions, the system releases or absorbs a fixed amount of energy while its temperature 
stays constant as heat is removed or added to the system. In this period, two different phases 
with similar temperatures and pressures but with two different internal energy densities co-exist, 
only their volume fraction changes, until one of them reaches 100 % and the other disappears. 
This phase coexistence region is frequently and incorrectly referred to as "mixed phase" although 
correctly this does not correspond to a single (but mixed) phase, instead it corresponds to a mixture 
of two, locally different, but co-existing phases. At very high baryon densities, the rehadronization 
phase transition of a quark-gluon plasma is believed to be such a first order phase transition, during 
which a specific amount of latent heat is released. As the energy density of the matter can change at 
a constant pressure, the speed of sound typically vanishes during such a first order phase transition. 
At the critical point {TcPLc) the latent heat disappears, and the order of the phase transition changes 
from a first order to a second order one. In the vicinity of the critical point, the difference between 
the gaseous and the liquid phase becomes almost non-existent. Fluctuations appear on all possible 
wavelengths and are signaled by critical opalescence, or the decreased attenuation length at all 
possible wavelengths. 

As the characteristic sizes disappear from the physical description of second order phase tran- 
sitions at the critical point, power-law distributions emerge and characterize the physical quantities 
and the exponents of these power-laws, called the critical exponents become the relevant, physically 
important observables. Frequently, a reduced temperature is introduced: 



c 



Tc 



(3.1) 



Along the line of phase transitions, a small change in temperature corresponds to a small change in 
baryochemical potential, Ar oc A/x^, so the power law exponents will be similar if the fluctuations 
are studied in terms of deviations from the critical temperature, or if t is re-defined as dimensionless 
deviation from the critical value of the baryochemical potential: t' = ^"^j^"'' ■ 

To characterize the phases, an order parameter |0| is introduced, which is zero before and 
non-zero after the phase transition, conventionally so that |(^>| > corresponds to T < 7^. In case of 
QCD phase transitions, a theoretical order parameter can be the expectation value of the quark con- 
densate, 1 1 = {qq) . Experimentally one can, for example, define an order parameter, that measures 
the violation of the quark number scaling of elliptic flow of identified kaons (vf ) and identified pro- 
tons (V2) , and use it as a phenomenological order parameter: |0'| = / djJL{w) [^vf (y) — ^V2(3^)] . 
For the definition of the universal scaling variable w and that of the elliptic flow V2 see ref. [p^], 
but note that at mid-rapidity w oc {mt — m)lT^, where mt — m = KEt is the transverse kinetic en- 
ergy and r* stands for the slope parameter of the single particle spectra. The limits of integra- 
tion are given to the upper kinematic limit of the applicability of a hydrodynamical picture, typi- 
cally pt ~ 1.5 — 2 GeV. The measure d}JL{w) may be chosen appropriately, for the simplest choice 
I recommend to try d}JL{w) = dw, but of course other measures like d}jL{\;w) = exp(— w)t/w or 
<ijU(A;w) = exp(— Aw)(iw etc are also possible. It is physically important to remove data points 
that are not from the hydrodynamic kinematic range, i.e. close to the target or projectile rapidity or 
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have large pt >2 GeV. The integral above has a vanishing argument if the elliptic flow corresponds 
to a perfect fluid of quarks, othervise it is positive, hence it can be used as an order parameter, chos- 
ing any form of the above measures. Of course similar order parameters can be defined for other 
particle pairs. Typically a baryon and a meson should be paired, to utilize their different number 
of valence quarks. I recommend to leave out pions, as their light mass (would-be Goldstone boson 
nature) and the unknown confinement of two heavy valence quarks into a single light pion makes 
their elliptic flow connection to the quark elliptic flow rather subtle and theoretically challenging. 

The new order parameter defined above is vanishing in a fluid of quarks and anti-quarks (better 
known as a strongly interactive quark-gluon plasma phase), while it is expected to become non- 
vanishing in a hadronic phase, i.e. when the quark number scaling of the elliptic flow is broken. 

We introduce below the six characteristic critical exponents of the second order phase transi- 
tions, a, 15, 7, 5, Tj and v following the lines of ref. |28]. 

The heat capacity at constant volume V or pressure p is defined as C = T ^\ where the deriva- 
tive is taken at constant values of V or p, respectively. Refs. [29, suggested that event-by-event 
fluctuations of the temperature T be used to measure the heat capacity at constant volume as 

1 (Ar)2 

7^ = H^. (3.2) 

A caveat is if the parameter T in this equation should be obtained from fits to the single parti- 
cle spectra (Ttm, the kinetic freeze-out temperature) or from fits to particle abundances {Tchem the 
chemical freeze-out temperature). Given the expectation that T in the above equation refers to tem- 
peratures in the vicinity of the critical point, I recommend to use temperature measures that yield 
values closer to the line of the phase transitions, e.g. T = Tchem, because Tchem > 71m- Note however 
that both temperatures may have a spatial distribution even in a single event, so this measure can 
be considered only as an approximate one. In a second order phase transition, the critical exponent 
a measures the power-law behavior of this heat capacity near the critical point: 

C(0=^ir", (3-3) 

where a may have both negative and positive values, as for example in the case of onset of super- 
fluidity in ^He, or for the case of the 3d Ising model, respectively. If the baryochemical potential 
tends to vanishing values as the critical point is approached, /ig = a{T — Tc)^ for some positive 
constants a and b, the specific heat capacity of baryons obeys the above law with a = 2, as shown 
in ref. [0]. 

Near but below the critical point, the temperature dependence of the order parameter is given 
by the critical exponent j8 : 

\<^{t)\c^\t\P , for t <Q. (3.4) 
The correlation function of the order parameter is defined as 

G(r,O = (0(r,O0(O,O)- (0(0,0)', (3-5) 
where r is a measure of relative coordinates. At large distances and for f / 0, 

G(r,0^^exp(-|r|/C), (3.6) 
\r\ 
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where A=A{t) and the correlation length ^ = ^{t) are both temperature dependent, but both are 
independent of relative separations. The susceptibility exponent 7 is defined as 

I (PrG{r,t) oc \t\-'^. (3.7) 

The exponent 5 in the 3d Ising model is related to a small external field H that breaks the 0(4) 
symmetry of the model explicitly. In terms of QCD, such a symmetry breaking term is related 
to the small but non-vanishing value of the masses of light quarks, H oc = niu = nid- In the 
presence of an external field, the order parameter is not vanishing at the critical temperature, but 
scales with the magnitude of the symmetry breaking terms as 

= 0,//^0))oc (3_8) 
Two very important critical exponents are V, the exponent of the correlation length ^, defined as 

c-ir', (3.9) 

and T], the exponent of the Fourier- transformed correlation function in the critical point {t = 0), 
defined as 

G(/t^O)oc (3.10) 

Out of these six critical exponents only two are independent. The remaining four are related by the 
theory of second order phase transitions as 

a = 2-dv, 

15 = X(j_2 + r7), 

7 = (2-r])v, 

<3 — 2 + T] 



where d stands for the number of dimensions, J = 3 for heavy ion collisions. See ref. [ ]28p and 
references therein for more details on how to obtain these relationships from the assumption that the 
free energy at the critical point is described by scale invariant, generalized homogeneous functions. 

It is quite remarkable that second order phase transitions of apparently very different systems 
are often characterized by the same set of critical exponents. This phenomenon is called universal- 
ity. For example, critical exponents at the critical point of the liquid-gas phase transition have been 
found to be independent of the chemical composition of the fluid. Furthermore, they are exactly 
the same as the critical exponents of the ferromagnetic phase transition in uniaxial magnets, which 
is the same as that of the 3 dimensional Ising model. Such systems are then referred to as being 
in the same universality class. The renormalization group theory of second order phase transitions 
states that the thermodynamic properties of a system near a critical point depend only on a small 
number of features such as the dimensionality of the system and the symmetry of the interactions, 
but are insensitive to the underlying microscopic properties. That is why it will be so interesting to 
test, if the perfect liquid of quarks, created in Au-i-Au collisions at RHIC, features a critical point 
or not. It is the universality class of QCD that characterizes matter at this point. This universality 
class is the key observable that experiments will have to determine, if indeed such a critical point 
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of QCD is found. Rajagopal and Wilczek predicted already in 1992 [28] that the universality class 
of QCD will be that of the 3d Ising model, where the critical exponents are measured as 



^7 
V 



0.03 ±0.01, 
0.73 ±0.02, 



hence the remaining exponents are given as 



(3.12) 



a 

7 
5 



-0.19 ±0.06, 
0.38 ±0.01, 
1.44 ±0.04, 
4.82 ±0.05. 



(3.13) 



We have conjectured that the violent environment of a heavy ion collision generates random 
fields that interact with the expanding medium and hence change the universality class of the 3d 



Ising model to that of the random field 3d Ising model |33]. In particular, due to this expected noise, 



the extremely low correlation exponent rj increases, and the critical exponents are given [ ]32[ ] as 



Tj = 0.50 ±0.05, 
V = 1.1 ±0.2. 



(3.14) 



hence the remaining exponents are given as 



a = -1.3±0.6, 
j8 = 0.6±0.1, 
7 = 2.2 ±0.4, 
5 = 4.7 ±0.3. 



(3.15) 



Note that the changes are most significant in the exponents 77 and v and these are the exponents 
that are most directly measurable, as will be discussed below. 



4. Soft correlations for measuring critical exponents of the correlation function 

In this subsection, we discuss how to measure the critical exponent of the correlation function 
77, based on a shape analysis of the Bose-Einstein or HBT correlation functions. 

Bose-Einstein correlation functions in high energy physics for Levy stable sources were intro- 



duced in refs. [34, B9p. These papers and ref. [33| are recited in this subsection, to point out how 



to measure the critical exponent of the correlation function in the vicinity of a critical point. 

Bialas considered a distribution of Gaussian coordinate space distributions, where the Gaus- 
sian radius parameter has a power-law distribution [^], and related the Bose-Einstein correlations 
of such sources to scale-invariant multiplicity fluctuations called intermittency. Brax and Peschan- 
ski introduced Levy distributions to multiparticle production in high energy physics [pj]], but in 
momentum space. They suggested to use the measured values of the Levy index of stability to sig- 



nal quark gluon plasma production in heavy ion physics. Here we follow ref. [ 33 1 in reconnecting 
these apparently different topics. 
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4.1 Bose-Einstein correlations for Levy stable source distributions 

The Bose-Einstein correlation function of identical boson pairs can be defined with the help 
of the corresponding two-particle and single-particle invariant momentum distributions: 

C2(ki,k2) = 4#4^. (4.1) 
M(ki)M(k2) 

If the long-range correlations can be either neglected or be corrected for, and the short-range cor- 
relations are dominated by Bose-Einstein correlations, this two-particle Bose-Einstein correlation 
function is related to the Fourier-transformed source distribution. 

For clarity, let us consider a one-dimensional and static source, where the space-time distribu- 
tion f{x) and the momentum space distribution g{k) are factorized: 

S{x,k)=f{x)g{k). (4.2) 

The corresponding Bose-Einstein correlation function is 

C2{kuk2) = l + \f{q)\\ (4.3) 

where the Fourier transformed source density (often referred to as the characteristic function) and 
the relative momentum are defined as 

f{q) = j ^ exp(/^x)/(x), q = ki-k2. (4.4) 

During a second order phase transition, where fluctuations appear on all possible scales with a 
power-law tailed distribution, the final position of a particle is given by a large number of position 
shifts, hence f{x), the distribution of the final position x is obtained as 

x = fix) = [ IVI^,dxin"j^Jj{xj)5{x- £x,). (4.5) 

!=1 k=l 

A system at a second order phase transition becomes invariant under a renormahzation group trans- 
formation. If we consider particle production near to such a scale invariant point, we may expect 
that adding an (n -|- l)th random shift does not change (after rescaling) the overall source distribu- 
tion. Thus at a second order phase transition we expect f{x) to be stable for convolution. 

In probability theory. Central Limit Theorems state that under certain conditions, the distribu- 
tion of the sum of large number of random variables converges (for n ^ oo) to a Umit distribution. 
If the elementary steps obey a "normal" elementary distribution that has a finite mean and variance, 
then the limit distribution of their sum can be proven to be a Gaussian. However, in the vicinity 
of a second order phase transition point, fluctuations appear on all scales and the variance of the 
elementary process diverges. In this case, the probability distribution will deviate from a Gaussian, 
in particular it develops power-law tails, but it still may be stable for convolutions. 

Stable distributions are limit distributions that occur in Generalized Central Limit theorems. 
The study of such distributions was started by the French mathematician Paul Levy in the 1920's. 
Stable distributions can be given simply in terms of their characteristic functions, as the Fourier 
transform of a convolution is a product of the Fourier-transforms, 

fiq)=flMq)- (4.6) 
(=1 
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C(q) BEC/HBT for Levy 




Figure 3: Bose-Einstein or HBT correlation functions for Levy stable source distributions for various indices 
of stability ■ At a second order QCD phase transition, the index of stability becomes one of the critical 
exponents of the phase transition, namely the exponent of the correlation function, Tj = Uuvy Note that on 
this plot, a stands for the Levy index of stability, au^-y, and it is not to be confused with the critical exponent 
of the heat capacity, which is denoted by the same symbol in the body of the manuscript. From Ref. [^. 



The characteristic function of univariate and symmetric stable distributions is 

/(^) =exp(/^5- 1 7^1 



(4.7) 



where the support of the density function f{x) is (— oo^oo). Deep mathematical results imply that 
the index of stability, au„y, satisfies the inequality < a^^,, ^ 2, so that the source (coordinate) 
distribution be always positive. These Levy distributions are indeed stable under convolutions, in 
the sense of the following relations: 



fi[q) = exp {iqSi - ImI"""") > ni;(9) = exp {iq5 - \yq\''^^^) , 



!=1 



(=1 

n 

5 = 1^5, 

i=\ 



(4.8) 
(4.9) 



Thus the Bose-Einstein correlation functions for uni-variate, symmetric stable distributions 
(after a core-halo correction, and a re-scaling) read as 



C{q\au„) = 1+Aexp(-|^/?| 



(4.10) 



See refs. [34] and [39] for further examples and details. 



4.2 Bose-Einstein correlations at a second order QCD phase transition 

At the CEP, the second order phase transition is characterized by the fixed point of the renor- 
malization group transformations. In a quark-gluon plasma, the vacuum expectation value of 
the quark condensate = l^q) vanishes, while in the hadronic phase, this vacuum expectation 
value becomes non-zero. The correlation function of the order parameter is defined as p(r) = 
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(0(r)0(O)) — (0)^ and measures the spatial correlation between pions. At the CEP, this correlation 
function decays as 

p(r)ocr-(^-2+Ti)^ (4_;^j) 

a power-law. The parameter T] is called as the exponent of the correlation function. 

For Levy stable sources, corresponding to an anomalous diffusion with large fluctuations in 
coordinate space, the correlation between the initial and actual positions decays also as a power- 
law, where the exponent is given by the Levy index of stability a^^^, as 

p(r)ocr"(i+«^-.v). (4.12) 

As we are considering a QCD phase transition in a = 3 three-dimensional coordinate space we 
find that the correlation exponent equals to the Levy index of stability, T] = auvy Fig. || indi- 
cates that the change in the shape of the correlation function is rather significant, if a^^,,, decreases 
from its Gaussian value of 2 to 0.5, its characteristic value conjectured from the random field 3d 
Ising model. Fig. ^ illustrates how this shape parameter of the Bose-Einstein correlation function 
may depend on the reduced temperature near the critical point. Hence the shape parameter of the 
two-particle Bose-Einstein or HBT correlation functions can be used to determine if the pions are 
emitted from the neighborhood of the critical end point of the QCD phase diagram [p3|]. Unfortu- 
nately, the bombarding energy or the centrality of this shape parameter has not yet been determined 
in heavy ion collisions. 

It is recommended that such a sharp drop of the index of stability a^yy is searched first as a 
function of the beam energy, the size of the colliding nuclei and the centrality of events, by mea- 
suring the excitation function of the Levy index of stability in the neighborhood of the maximum 
of opacity K. Close to the critical point, where such a drop in the Levy index of stability is already 
observed, the collected events could be further separated to several sub-classes based on the event- 
by-event values of Tchem, the chemical freeze-out temperature, fitted to the chemical composition 
of the given event. This way, the dependence on T = (T — r( )/rc of observables Uke the critical 
exponent of the correlation function or that of the correlation length can be measured. 

4.3 Negative binomial multiplicity distributions and the exponent of the correlation length 

In this subsection we discuss how to measure the critical exponent of the correlation length T] . 
An indirect method to measure a quantity proportional to the correlation length has already been 



published by PHENIX in ref. [|42|]. Observing that the Bjorken energy density estimate, hence the 
expected initial temperature, scales as a monotonous function of the number of participants, that pa- 
per studied the centrality dependence of the multiplicity distributions in y^S^ = 200 GeV Au-i-Au 
collisions at RHIC. In various intervals of pseudo-rapidity T] = 0.5 log (^|^p^)> the multiplicity 
distribution was fitted by a negative binomial distribution (NBD) 

that has two fit parameters, the mean multiplicity {n) and the width parameter k. This parameter is 
related to the second scaled factorial moment of the multiplicity distribution, 

1 {n{n-l)) 



k {nY 



1. (4.14) 
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Levy index of stability 




Figure 4: To locate and characterize the QCD critical point (the endpoint of the line of the assumed first 
order QCD transitions at high baryon densities) the Levy index of stability a = Uuvy is plotted as a function 
of reduced temperature, the dimensionless variable t. From Ref. 



Data were tabulated in one dimension, as a function of the pseudo-rapidity. For sufficiently large 
pseudo-rapidity intervals, where the correlation length is smaller than the size of the window, <^ << 
Sr], the width parameter is proportional to the correlation length: 



k or\ 



(4.15) 



and using this property the product al^ has been determined. Most interestingly, an indication 
was found for a non-monotonous behavior of the correlation length at centralities corresponding to 



about 90 participating nucleons [42], in the 35% - 45% centrality class 



However, other observables including opacities of Table ^ indicate no signals for a non- 
monotonic behavior in the vicinity of this centrality class. Clearly, more detailed studies are nec- 



essary, the opacity estimates of Table ^Jj for example depend on the assumption that the Bose- 
Einstein correlation function can be fitted with a Gaussian form, hence a Gaussian Rhbt is a good 
measure of the length-scale over which the attenuation happened. However, such an assumption is 
clearly not adequate when looking for a critical point. A possible decrease of the shape parameter 
aLn^y of the correlation function at this centralities should be investigated together with the question 
if a Levy fit to the Bose-Einstein correlation function would yield an improved description of Bose- 
Einstein correlation data and if the extracted scale parameters /?i^„,, would be significantly different 
from the Gaussian radius parameters. 

Similarly to the suggested study of the exponent, a large sample of events, collected close to 
the interesting centrality class of 25- 35 % , need to be separated to several sub-classes based on the 
event-by-event fluctuating values of Tchem, the chemical freeze-out temperature, given by the hadro- 
chemical composition of each event separately. This way, the dependence on T = {T — Tc)/Tc of 
the correlation length and that of the opacity could be measured, to see if indeed in ^/sNN = 200 
GeV Au-i-Au collisions a critical phenomenon at 25-35 % centralities is approached, or not. 
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5. Characterizing a second order phase transition in QCD: a strategy 

The following straightforward strategy is suggested to locate and characterize experimentally 
the properties of a second order QCD phase transition: 

a) Locate experimentally the critical end point (CEP) of the line of 1st order phase transitions 
by searching for critical opalescence, i.e. looking for the maximum of opacity, or the minimum of 
attenuation length, as a function of centraUty, bombarding energy and size of the colUding nuclei 
in a broad, high transverse momentum range. 

b) Determine the chemical freeze-out temperature T on an event-by-event basis for an event 
selection that corresponds to the maximum of opacity. Check that various observables (heat ca- 
pacity, correlation length) indeed indicate a power-law behavior as a function of t = \T — Td/Tc 
simultaneously for the event class selected in point a). 

c) Determine the critical exponents for at least two observables. Use theoretical relations to 
determine the remaining critical exponents, so that all the six critical exponents be given. It is 
recommended to try to measure the critical exponent of the correlation function rj as it is defined in 
the momentum-space, and the power-law shape shows up as a function of the relative momentum 
and not as a function of t, hence this exponent is easy to measure. The correlation length can 
be determined from the multiplicity distribution using the negative binomial distribution, and the 
exponent of the correlation length V can be obtained from event-by-event determination of both the 
correlation length and the chemical freeze-out temperature T. 

d) Determine the universality class of QCD from the two measured critical exponents. 

e) Cross-check the result by measuring additional critical exponents and cross-check these 
measurements with constraints obtained from the theory of second order phase transitions. For 
example, determining the exponent of the heat capacity a seems to be straight-forward, as the heat 
capacity is measured by the event-by-event fluctuations of the chemical freeze-out temperature. 

6. Conclusions 

The results summarized here are relevant for the low energy program of the RHIC accelerator, 
to second or third generation experiments in the heavy ion programme at CERN SPS as well as for 
the cold and dense hadron matter physics planned at the forthcoming FAIR facility. The methods 
outlined can provide means to locate, if it exists, the critical end point of the line of first order phase 
transitions in QCD. Based on an analogy with experiments in physical chemistry and statistical 
physics, a novel method to determine opacity was proposed and critical opalescence was suggested 
as a smoking gun experimental signature for the vicinity of the critical point. By measuring the 
critical exponents of the correlation function and the correlation length with the outUned methods, 
the universaUty class of the QCD critical point can be determined. 
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